Only n-Qubit Greenberger-Horne-Zeilinger States are Undetermined by their 

Reduced Density Matrices 
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The generalized n-qubit Greenberger-Horne-Zeilinger (GHZ) states and their local unitary equiv- 
alents are the only states of n qubits that are not uniquely determined among pure states by their 
reduced density matrices of n — 1 qubits. Thus, among pure states, the generalized GHZ states are 
the only ones containing information at the n-party level. We point out a connection between local 
unitary stabilizer subgroups and the property of being determined by reduced density matrices. 

PACS numbers: 03.67.Mn,03.65.Ta,03.65.Ud 



Quantifying and characterizing multi-party quantum 
entanglement is a fundamental problem in the field of 
quantum information. Roughly speaking, one expects 
the states that are "most entangled" to be the most valu- 
able resources for carrying out quantum information pro- 
cessing tasks such as quantum communication and quan- 
tum teleportation, and to give the most striking philo- 
sophical implications in terms of the rejection of local 
hidden variable theories [ij. 

Although no single definition of "most entangled" 
seems possible, since we know that multi-party entangle- 
ment occurs in many types that admit at best a partial 
order 0] , it is still worthwhile to consider properties that 
carry some of the spirit of "most entangled." 

One such property is the failure of a state to be de- 
termined by its reduced density matrices. As reduced 
density matrices contain correlation information pertain- 
ing to fewer than the full number of parties in the sys- 
tem, states exhibiting entanglement involving all parties 
must possess information beyond that contained in their 
reduced density matrices. Linden, Popescu, and Woot- 
ters put forward this suggestion in [sl, [j| and proved the 
surprising result that almost all 7i-party pure states are 
determined by their reduced density matrices. In other 
words, the set of n-party pure states undetermined by 
their reduced density matrices is a set of measure zero. 
In Q , Diosi gave a constructive method that succeeds in 
almost all cases for determining a 3-qubit pure state from 
its reduced density matrices. Nevertheless, the question 
of precisely which states are determined by their reduced 
density matrices remained open. 

In this Letter, we show that the only n-qubit states 
that are undetermined among pure states by their 
reduced density matrices arc the generalized n-qubit 
Greenberger-Hornc-Zcilingcr (GHZ) states, 

a|00---0)+/3|ll---l), a,/3^0 

and their local unitary (LU) equivalents. This means 
that, among pure states, the generalized GHZ states are 
the only ones containing information at the n-party level. 
For the case n = 3, this result was reported previously in 



[3|. Part of our argument employs the methods of ,5| in 
an essential way. 

Let Dn be the set of n-qubit density matrices. If p G 
Dn is an n-qubit density matrix, and j € {1, . . . , n} is a 
qubit label, we may form an (n— l)-qubit reduced density 
matrix = trj p by taking the partial trace of p over 
qubit j. Let 



PTr : D,„ 



be the map p i-^ ■ • ■ i P{n)) that associates to p its 
n-tuple of (n — l)-qubit reduced density matrices. The 
map PTr is neither injective (one-to-one) nor surjective 
(onto). Its failure to be surjective means that there are 
n-tuples of (n— l)-qubit density matrices that cannot be 
produced from any n-qubit density matrix by the partial 
trace. The question of whether a collection of (n — 1)- 
qubit reduced density matrices could have come from an 
71-qubit density matrix by the partial trace is the subject 
of recent and ongoing investigations @, 0]. The failure 
of PTr to be injective means that multiple n-qubit states 
can have the same reduced density matrices. States pi ^ 
P2 with PTr(pi) = PTr(p2) require more information for 
their determination than is contained in their (n — 1)- 
qubit reduced density matrices. 
Let 



Pr, 



{peD.n\p^ 



P} 



be the set of pure n-qubit states. If we are interested 
primarily in pure states, we can restrict the partial trace 
map to pure n-qubit states. 



ptr = PTr|p„ -.Pr. 



D" 



Given a pure state tp with {ip\ g P„, the set 
ptr^^(ptr(?/;)) contains all pure states with the same re- 
duced density matrices as ip. (We abbreviate ptr(?/') — 
ptr(l^) 

We define a state "0 to be determined among pure states 
if ptr~^(ptr(-0)) contains only and undetermined 

among pure states if ptr~^(ptr (?/;)) contains more than 
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one state. Similarly, we define a state p € D„ to be de- 
termined among arbitrary states if PTr~^(PTr(p)) con- 
tains only p, and undetermined among arbitrary states if 
PTr~^(PTr(p)) contains more than one state. 

The surprising result of Linden and Wootters [1] is 
that almost all n-qubit pure states are determined among 
arbitrary states. 

Nevertheless, there are pure states that are undeter- 
mined among pure states (and consequently undeter- 
mined among arbitrary states). For example, consider 
the onc-paramctcr family of n-qubit states 

l^)^_L|oo...o) + -||ii...i), 

where 77 is a complex number with magnitude one. If i]i 7^ 
772, then 1 771) and \r]2) are different states with different 
density matrices I771) (r/ij ^ I772) (772 1, yet they share the 
same reduced density matrices, that is, ptr(|7yi) ~ 
Ptr(|772) (?72|)- 

We sec that almost all pure 7i-qubit states are deter- 
mined among pure states, yet 7T,-qubit GHZ states are 
undetermined among pure states. The question then be- 
comes, precisely which states "0 are undetermined among 
pure states? 

Main Result. An 7i-qubit state 7/; is undetermined 
among pure states if and only if -0 is LU-equivalent to 
a generalized n-qubit GHZ state. 

Proof. Let ij: be an n-qubit pure state. 

Suppose that ip is LU-equivalent to a generalized n- 
qubit GHZ state, so we have 

[/|7A) ==a|00---0)+/3|ll---l), 

where [/ is a local unitary transformation. Define 
f/|V'') = a|00---0) - Then It/;) (t/^I ^ 

(?/>'|, since a/3 7^ 0, but ptr(7/;) = ptr(7/''). Hence, 
Ip is undetermined among pure states. 

Conversely, suppose that 1-0) is undetermined among 
pure states. Then there is an n-qubit state vector ^ 
e*" It/i) that has the same reduced density matrices as 

Claim: If [t/;) and have the same reduced density 
matrices, then for each qubit j G {!,..., n}, there is 
a one-qubit local unitary transformation Lj such that 
W)=h\^P). 

To prove this, let j e {1, . . . , n} be a qubit label. Let 
Pj denote the one-qubit reduced density matrix of It/;) for 
qubit ). We write pj as a spectral decomposition. 



The (n — l)-qubit reduced density matrix 

P{3) = ^^J IV') (V'l , 

obtained by taking the partial trace of \^) (t/;| over qubit 
J, has the same nonzero eigenvalues as pj. 



If 7^ pj, then the (n— l)-qubit eigenvectors |0; {j)) and 
|1; ())) arc unique up to a phase. If ^ -P]j then the 
eigenvectors of p^) with eigenvalues^ ~ P] ~ ^1'^ consti- 
tute a two-dimensional subspacc of the ^-dimensional 
vector space of (n — l)-qubit vectors, and any orthonor- 
mal pair of vectors in this subspacc may be chosen as a 
basis. 

We choose the one-qubit orthonormal basis \ij) and 
the (n — l)-qubit orthonormal basis \ij] (j)) so that |t/') 
can be written 

|tA> = |0) ®, |0; (j)) + yp] |i) ®, |i; (j)> , 

where (^j is the tensor product that inserts a one-qubit 
ket just before the jth factor in the (n — l)-qubit ket 

K.;(j)>- 

Now can be regarded as the state of a bipartite 
system composed of qubit j and all qubits but j , and has 
a Schmidt decomposition with respect to those subsys- 
tems. 



m ®,|0';(j)> 



|i'> ®.|i';(j)>, 



where |0'),|1') are orthonormal one-qubit vectors and 
jo'; (j)) , |1'; (j)) are orthonormal (71 — l)-qubit vectors. 
Taking the partial trace over qubit j, we have 



tr, IV') (V-'l 



^"lO';^-)) (0';(j)| 



-qj|l';(j)) (l';(j)|. 



Since this must be equal to it must have the eigen- 
values of /9(j), q'j = p'j and ~ p)- We consider two 
cases, depending on whether p^j) has distinct eigenvalues 
or not. Let us treat first the case of distinct eigenvalues, 
p° 7^ p]- t^^i^ case, the eigenvector |0'; [j)) can be 
off by at most a phase from the eigenvector |0; (j)), and 
similarly for |1'; {])) . The same argument applied to the 
one-qubit reduced density matrix pj shows that |0') can 
be off by at most a phase from |0), and similarly for |1'). 
In this case, then, we can write 



Pj = 



for some orthonormal basis \ij), where p'j and arc the 
eigenvalues of pj . If 7^ pj , then the orthonormal basis 
\ij) is uniquely determined up to a phase. If p" = pj, 
then any one-qubit orthonormal basis can be used. 



/pO(L,|0))®, |0; {3)) + Jp]{L,\l))> 



|i;W>, 

(1) 

with a 2 X 2 diagonal unitary matrix. 

Let us treat next the case of repeated eigenvalues, 
p° = pj. In this case, the eigenvectors |0'; (j)) and 
|1'; {j)) must merely span the same two-dimensional com- 
plex space that is spanned by |0; (j)) and |1; ())). In this 
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case, the primed eigenvectors must be related to the un- 
primed eigenvectors by a two-dimensional unitary trans- 
formation, 

\Q';{j))=uoo\0;U))+uoi\l;{j)) 



with 



"00 "01 
"10 "11 



G U{2). 



The same argument applied to the one-qubit reduced 
density matrix pj shows that there must be some 2x2 
unitary matrix vim with 

|0') ="oo|0)-t-«oi|l) 
|1'> ="10 |0)+ Willi). 

In this case, then, we can write equation ([T]) with Lj the 
2x2 unitary matrix equal to the product of the transpose 
of Vim with uirn- (Wc havc abused notation by using the 
symbol Lj to represent both the 2x2 unitary matrix, 
and also the local unitary transformation on n-qubit state 
vectors 



/ (g) • • • (8) / (8) «) / ( 



with the 2x2 matrix Lj in the jth slot of this tensor 
product, and onc-qubit (2 x 2) identity operators in all 
other slots.) This completes the proof of the Claim. 
For each pair of qubit labels j, k, we have 

1^) = L^'L, W . 

Next, spectrally decompose each Lj with unitary matri- 
ces Uj so that 

Dj = U,L,UJ' 
are diagonal. Wc have 

D^^DjUi •••[/„ IV^) = C/i • • • UnU^'D^'UkU-'DjUj |V) 

Using the multi-index / = (^1^2 ■ • • in), where each ij is 
zero or one, and the basis 



|/) 



expand 



\^l^2 = \n) <X) \i2) 



Ui---Un\yj)=J2ci\I) 



where e'^^ 7^ e , since \ip') is a different state from 
lip). Now we havc 

CI = cicxp{i[aj -ak + f3j - {-1)'" l3k]} 

for all multi-indices / and all j. A; e {1, . . . ,n}. We see 
that for each multi-index /, either c/ = or 

exp{i[aj - ak + (-l)'^/3j - {-1)'" Pk]} = 1 

for all j, fc e {l,...,7i}. Let J be a multi-index with 
cj 7^ 0. If / is any multi-index that agrees with J in at 
least one entry (say the jth qubit entry), and disagrees 
with J in at least one entry (say the fcth qubit entry), 
then c/ = 0. We conclude that cj and Cj, where J is the 
multi-index consisting of the complements of each of the 
n bits in multi-index J, are the only nonzero coefficients 
in Ui ■ ■ ■ Un \tp)- Consequently, \tp) is LU-equivalent to a 
generalized n-qubit GHZ state. □ 

Much of this argument carries over to the more general 
situation of a system of n parties in which party j has 
dimension dj (party j is a qubit if and only if dj = 2). In 
particular, the Claim carries over. Suppose that \ip) and 
lip') are states of a system of n parties in which party j 
has dimension dj. If \tfj) and have the same reduced 
density matrices, then for each party j S {1, . . . , n}, there 
is a one-party local unitary transformation Lj such that 
\tl)') = Lj The second part of the argument is com- 
plicated by the possibility of repeated eigenvalues in the 
transformations Dj. We leave this as a question for fu- 
ture work. 

It is worthwhile to point out the significance of stabi- 
lizer subgroups of the local unitary group in this work. 
The local unitary group for n-qubit density matrices is 
the group G = S'J7(2)", consisting of a special unitary 
transformation on each qubit. Each (pure or mixed) state 
p has a stabilizer subgroup Ip consisting of elements of G 
that leave p fixed under the action gpg~^ for g € G. We 
have seen that a state that is undetermined among pure 
states has a special type of enlarged stabilizer subgroup. 

There are two alternative formulations of the main re- 
sult that may suggest promising avenues for the classi- 
fication of entanglement types. One can precisely char- 
acterize the pure n-qubit states that are undetermined 
among pure states in terms of the stabilizer subalgebra 
of the state, that is the Lie algebra of the stabilizer sub- 
group of the state. It is shown in Q that the generalized 
n-qubit GHZ state has (for n > 3) stabilizer subalgebra 




and write 



e«ft 
e-'P^ 



where Zj is the Pauli matrix tr^ applied to qubit j. States 
of n qubits undetermined among pure states are precisely 
those which are LU-equivalent to states with this subal- 
gebra. 
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A second alternative formulation of the main result is 
in terms of the dimension of the stabilizer subgroup. For 
n = 3 and n> 5, an n-qubit state is undetermined among 
pure states if and only if it is not a product state and its 
stabilizer subgroup has dimension n — 1 0] . 

We have shown that all n-qubit states other than gen- 
eralized n-qubit GHZ states and their LU-equivalents are 
completely determined by their reduced density matrices. 
Is it necessary to specify all of the reduced density matri- 
ces? Which states are undetermined by specifying only 
n — 1 (rather than all n) of the {n — l)-qubit reduced 
density matrices? Is it a larger set than the generalized 
GHZs? The answer is yes, and stabilizers can help us 
understand this. For example, the state 

|X) = ^(10000) + 10001) + 11111)) 

is undetermined by its 3-qubit reduced density matrices 
obtained by taking the partial trace over qubit 1, the 
partial trace over qubit 2, and the partial trace over qubit 
3. It is not LU-equivalent to a generalized 4-qubit GHZ 
state (it has a different stabilizer subalgebra structure, 
and stabilizer subalgebra structure is an LU-invariant). 
Note that Zi Ix) = ^2 Ix) = ^3 Ix) 7^ e*" |x). The state 
Zi Ix) has the same 3-qubit reduced density matrices as 
Ix) when taking the partial trace over qubit 1, 2, or 3. 
These two states have a different 3-qubit reduced density 
matrix when taking the partial trace over qubit 4. 

A pure state's LU-equivalence class is often considered 
to contain all of the information about the entanglement 
of the state. An interesting question is: Are there n-qubit 
pure states with entanglement information that is not 
contained in their reduced density matrices? We might 
interpret this question as equivalent to the question: Are 
there n-qubit pure states for which the LU-equivalence 
class of the state is undetermined by its reduced density 



matrices? The answer to this question is no. Every n- 
qubit pure state can be determined (among pure states) 
up to a local unitary transformation by its reduced den- 
sity matrices. This can be seen directly from the Claim 
at the beginning of our proof. Two pure states that have 
the same reduced density matrices must be LU equiva- 
lent. 

We have not answered the question of which n-qubit 
pure states are undetermined among arbitrary states by 
their reduced density matrices. The set of n-qubit states 
undetermined among arbitrary states must contain the 
generalized GHZ states, bit it could be strictly larger 
than the set that is undetermined among pure states. 
This remains an open question. 
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